Graphene-based voltage-tunable coherent terahertz 

emitter 
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Abstract — A portion of the electromagnetic wave spectrum 
between ~ 0.1 and ~ 10 terahertz (THz) suffers from the 
lack of powerful, effective, easy-to-use and inexpensive emitters, 
detectors and mixers. We propose a multilayer graphene - boron- 
nitride heterostructure which is able to emit radiation in the 
frequency range ~ 0.1 — 30 THz with the power density up 
to ~ 0.5 W/cm^ at room temperature. The proposed device is 
extremely thin, light, flexible, almost invisible and may completely 
cover the needs of science and technology in the sources of 
terahertz radiation. 

Index Terms — Terahertz radiation, emitters, graphene, boron 
nitride. 

IF a fast electron moves in a periodic potential U{x) ^ 
sin(27rx/a:r) with the average velocity vq, its momentum 
Px, as well as velocity v^ = Px/m, oscillates in time with the 
frequency 

(1) 






Since electrons are charged particles, such a motion is accom- 
panied by an electromagnetic radiation with the frequency ([T]), 
Ref. fV\. This physical principle is used in backward- wave 
oscillators and free-electron lasers, where the potential U{x) 
is produced by periodic in space electric or magnetic fields. 
The frequency of radiation can be tuned in these devices by 
varying the accelerating voltage which determines the electron 
velocity vq. 

It was a long dream of scientists and engineers to create 
a compact solid-state emitter operating on the same phys- 
ical principle. For example, ||2l-(|5] such an idea could be 
realized in semiconductor structures with a two-dimensional 
(2D) electron gas: placing a metallic grating in the vicinity 
of the 2D conducting layer and driving electrons across the 
grating stripes it seemed to be possible to force all electrons 
to emit electromagnetic waves at the frequency ([T]i. However, 
instead of the strong coherent emission at the frequency ([T} 
a weak thermal radiation at the frequency of two-dimensional 
plasmons 



fp 



Uge^ 



m'eax 



(2) 



was observed |l2l-||5|; here e, to* and Ug are the charge, the 
effective mass and the surface density of 2D electrons and e 
the dielectric constant of surrounding medium. 

The reason of this failure was explained in Ref. 16|. It 
was shown that the single-particle formula ^ is valid only 
^t fp <C vq/gx, i.e. if the density of electrons is low or 
the drift velocity vq is sufficiently high. This condition is 



easily satisfied in vacuum devices and free-electron lasers. In a 
dense solid-state plasma one should take into account electron- 
electron interaction effects. Then one finds [6| that the strong 
coherent radiation is observed at the frequency 



f= — -fp 
a-x 

only if the velocity vq exceeds a threshold value 



vo > vth — fpax 



Hse^Qx 



(3) 



(4) 



for details see ||6|. Otherwise, at vq -C fpUx, electrons emit at 
the plasma frequency (|2]i, just due to the heating of the system 
(the thermal radiation). 

The condition (|4|i is of crucial importance for successful 
device operation. It shows that the surface electron density 
Us should be low and the drift velocity should be large. 
This condition is very difficult, if possible at all, to fulfil in 
semiconductor structures. For example, in a GaAs quantum 
well with TO* = 0.067TOO, e = 12.8, n^ - 3 x 10" cm^^ 
and Ux ^ 1 /im the threshold velocity (JUl is about 10* cm/s 
which is more than 4 times larger than the Fermi velocity at 
the same density. 

In graphene, a recently discovered ||2]-|l9] only one-atom- 
thin, truly two-dimensional carbon material, the spectrum of 
electrons is linear. 



^±(p) = ±VF 



pI 



Pi 



(5) 



with a very large, as compared to semiconductors, Fermi 
velocity vp ~ 10* cm/s. This makes graphene an ideal 
candidate for the realization of the old idea to create a voltage 
tunable solid-state terahertz emitter, see illustration in Figure 

m 

However, an attempt to directly apply the above described 
ideas to graphene faces two difficulties. First, in graphene 
oscillations of the momentum do not lead to oscillations of 
the velocity and, hence, of the current, since the velocity 
v± = 9i?±(p)/c'p is not proportional to the momentum. 
Second, the condition vth < vo < vp restricts the required 
electron density by the values of order of Ug ~ 10^° cm~^. 
It is known, however, that due to internal inhomogeneities 
(electron-hole puddles |10|) the average density of electrons 
in a graphene sheet typically exceeds ~ 10^^ cm^^ even at 
the Dirac point. 

The both difficulties can be overcome by using, instead 
of continuous graphene layers, an array of narrow stripes of 



Electrons in equilibrium 
Electrons in driving electric field 





momentum 

Fig. 1. In order to accelerate a small amount of electrons up to velocities 
of order of vp ft: 10* cm/s one needs a lai'ge change of the momentum in 
the case of the parabolic energy dispersion (semiconductors, left) and a much 
smaller change of the momentum in the case of the linear energy dispersion 
(graphene, right). 



graphene. This leads to the following basic design of the 
emitter 111], Figure |2] The first (active) graphene layer lies 
on a substrate made out of a dielectric material, e.g. Si02 
or hexagonal boron nitride (/i-BN), Figures [2^. This layer 
consists of an array of stripes with the width Wy and the period 
ay and has two metallic contacts "source 1" and "drain 1", 
Figure |2};. Above the first graphene sheet lies a thin dielectric 
layer, made out of a few monolayers of h-BN. A second 
graphene layer. Figure |2}l, has the shape of a grating with 
the stripe width Wx and the period ax, oriented perpendicular 
to the stripes of the first layer and covers the whole structure. 
It has a metallic "gate" contact. Figure |2}l. The side view of 
the whole structure is shown in Figure |2j). 

Due to the finite width Wy of graphene stripes in the first 
layer, the pj^ -component of the momentum is quantised, py ~ 
TTTi/Wy, and a gap is opened up in the graphene spectrum. 



Figure [3^. The value of the gap, 

_TThVF 



2An = 2- 



Wy ' 



(6) 



(7) 



is controlled by the stripe width Wy and can be as large as 
~ 40 meV (~ 500 K) if Wy =i 0.1 ^lm. The Hnear density of 
electrons Ue and holes Uh per unit length then reads 



ne{n,T) 



gs9v ■\p /"^ dx_ 

^^y n^l-'O 1 + CXp^^ T/Ao 



nh{tJ.,T) = ne{~^i,T), 



(8) 



where /i is the chemical potential, T is the temperature and 
<7s = <7u = 2 are spin and valley degeneracies. As seen from 
Figure|3}', the total linear charge density ?i; = rie + nh is about 
~ YljWy at room temperature. At Wy ~ 0.1 /im this gives 
ni « 1.2 X 10® cm^^. If to choose the period of the structure 
in the first layer ay bigger than 1.2 /im, the average surface 
density n^ ~ ni/ay will be smaller than 10^° cm~^. Thus, 
by choosing a sufficiently large ratio ay/Wy one can always 
satisfy the threshold condition ^. This can also be seen from 




Fig. 2. (a) The overall view of the device design. The first graphene layer 
lies on a substrate (made out of h-BN or Si02) and is covered by a few- 
nanometer thin /i-BN dielectric layer. Two metallic contacts, "source 1" and 
"drain 1", are attached to the graphene layer 1 from the "west" and "east" 
sides. On top of the ft-BN layer lies the second graphene layer with a metallic 
contact ("gate") attached on the "north" (or "south") side, (b) The side view 
of the graphene based emitter (c) The design of the 1st graphene layer The 
central (operating) area of the layer is made in the form of a periodic array 
of narrow stripes, (d) The design of the 2nd graphene layer (grating), (e) 
An alternative design of the top (grating) graphene layer with two contacts 
"source 2" and "drain 2". (f) A single-layer device structure with a modulated 
width of graphene stripes. 



the expression 



Vth 
VF 



mWy ttx 



ehvi 



(9) 



which should be smaller than unity (Eq. (|9]l follows from dUl 
if to replace m* by IS.q/v'^). The first factor e^ /ehvp here 
is of order unity. The second factor niWy depends on the 
temperature. Figure |3j), and is about 10 at room and about 
1 at liquid nitrogen temperature. The last factor ax /ay can 
be made as small as desired by the corresponding choice of 
geometrical parameters of the structure. 

In the operation mode. Figure ^p, a large dc voltage Vsd is 
applied between the source and drain contacts of the first layer 
and a small dc (gate) voltage V12 - between the first and the 
second graphene layer The source-drain voltage causes a dc 
current in the first-layer stripes flowing in the x-direction. The 




i'fPa/Ao 



Ai/Ao 



Fig. 3. (a) The energy spectram in the graphene stripes. If Wy ~ 0.1 /im, 
the gap 2Ao is about 500 K. (b) The linear charge carrier density ri; = 
^e + fih as a function of the chemical potential fi/Ao at different values 
of the temperature: T/Aq = 0.01 (dotted, blue), 0.3 (dashed, red), 0.8 (dot- 
dashed, magenta) and 1.2 (solid, black). If Wy ~ 0.1 fim, the blue, red, 
magenta and black curves approximately correspond to 2.5, 75, 200 and 300 
K, respectively. 



corresponding drift velocity should lie in the window vth < 
vq < vp. If the drift velocity is about 6 x 10^ cm/s and the 
grating period a^ — 0.2 /.tm, the fundamental frequency of 
the emitter will be about 3 THz. Dependent on the parameters 
this frequency can be red- or blue-shifted. 

Let us estimate the power of the emitted radiation. Under the 
action of the dc source-drain voltage electrons move with the 
constant drift velocity wq- A. small gate voltage V12 applied 
between the first and the second graphene layers results 
in the periodic time-dependent modulation of the velocity, 
i^xit) = vo + viCos{2TTV{)t/ax)- The modulation amplitude 
vi is controlled by the gate voltage V12 and can be varied 
between and vq. Assume that we have a 70% -modulation of 
the drift velocity, i.e. wi ~ 4x 10^ cm/s. Then the amplitude of 
the ac electric current density is ji ~ eUgVi, the electric and 
magnetic fields of the emitted wave (calculated from Maxwell 
equations) are £'1 = iJi ~ 27rji/c and the intensity of the 
emitted radiation (the Poynting vector) Wrad = cEI/Att. At 



the electron density rig 
radiation power of 

Wrad ^ - 



10 cm this gives a very large 



-{erisVi) 



0.5 W/cm^ 



(10) 



Estimating the Joule's heating power as Whe 



JqEq 



Jo/'^min, where a,nin = 4:3 /h is the minimal conductivity 
of graphene lE), we get 



Wr. 



Whe 



Jq I ^min 



< 



nar, 



he' 



(11) 



i.e. the efficiency of the emitter is about 1%. If the structure 
lies on a Si02 or BN 1 mm-thick substrate and the opposite 
side of the substrate is maintained at room temperature, the 
increase of temperature of the structure at Wheat — 50 W/cm^ 
does not exceed 10-20 K, due to the very large surface-to- 
volume ratio in the two-dimensional graphene and to the high 
thermal conductivity of the substrate. 

The proposed device may emit radiation not only at the 
fundamental frequency ^ but also at its harmonics. In the 
above discussion we have assumed that the periodic potential 
U{x) produced by the voltage V12 between the two graphene 
layers has a simple sinusoidal form U{x) — Uosm{2TTx/ar^). 
This is a good approximation for the case when the distance 





Vpl/Of harmonic number 

Fig. 4. (a) The time dependence of the velocity v{t) calculated for C/o/^O = 
0.3, ax'f/vp = 1, eBo(ii/Ao = 1 and S = 0.1 (black solid curve) and 
5 = 10 (red dashed curve). The average value of v(t)/vp in the cases S = 
0.1 and 5 = 10 are 0.661 and 0.615, respectively, (b) T he c orresponding 
Fourier spectra. Inset: The model periodic potential U{x) H4\ at 5 = 0.1 
(black solid curve) and 5 = 10 (red dashed curve). 



D between the 2D electron layer and the grating is compa- 
rable with the grating period. In the proposed graphene-BN- 
graphene device the thickness of the dielectric BN layer can 
be as small as a few nanometers while the realistic value of 
the grating period is > 0.1 /im. The ratio D/ax is therefore 
much smaller than unity and the potential U{x) has a step- 
like form, see Inset in Figure |4]3. The higher spatial harmonics 
will then lead to higher frequency harmonics in the emission 
spectrum. To show this quantitatively, we solve the Newton 
equations of motion for an electron in the graphene stripe 



—7- = -eEo - 7Pa; - 



dp. 
dt 
dx 
It 



dU(x) 
dx 



VpPx 



V^ 



vIpI 



(12) 
(13) 



where £'0 = Vsd/ L^ is the dc electric field caused by the 
source-drain voltage, L^ is the distance between the source 
and the drain and 7 = 1/r is a phenomenological parameter 
describing the "friction force" due to the scattering. For the 
periodic potential U{x) we use a simple model expression 



U{x) = [/o 



tanh[S' sin(27ra;/a2;)] 



(14) 



tanhS" 

where [/q is the potential amplitude and the parameter S 
determines the steepness of U{x), Figure |4j5: if S* < 1, the 
periodic potential is smooth and close to a single-harmonic 
sinusoidal form; if 5 ^ 1, it has a step-like form and 
contains many spatial Fourier harmonics. Equations (ITZt - 
(fT3T l are written for an electron of the first subband n = 1; 
evidently, taking into account the higher electron and all hole 
subbands will not significantly change the results. In these 
equations we also ignore the plasma effects since we assume 
that the electron density is sufficiently low, the parameter 
vth/vp, Eq. (|9]l, is small and the threshold condition (|4]i is 
satisfied, as has been discussed above. There are two sources 
of nonlinearity in equations (fT2li-(fT3]i: the non-sinusoidal peri- 
odic potential U{x) and graphene-specific nonlinear velocity- 
momentum relation. 

Figure |4] shows results of the numerical solution of equa- 
tions (fT2]i-(fT3]l for parameters Ua/Ao = 0.3, eEoax/Ao = 1 
and a^j/vp = 1- If 2Ao ~ 500 K ~ 40 meV the first two 
numbers correspond to Cq — 6 meV and E^ ~ 1000 — 200 
V/cm for the grating period a^ ~ 0.2 — 1 /im respectively. 
The last number means that the mean free path is of order 



of the grating period, i.e. ~ 0.2 — 1 /im for the same 
values of a^ (notice that the parameter cj/7 ~ lot is not 
very crucial: the device will work both at lot > 1 and at 
UJT < 1). It is quite possible to realise all these requirements 
in experiments. Figures |4^ and |4j5 show the time dependence 
of the velocity v.j;{t) and its Fourier spectra respectively; 
the black (solid) and the red (dashed) curves correspond to 
the small {S = 0.1) and large {S — 10) values of the 
steepness parameter S. One sees that at S* <C 1 the current 
jx{t) = eUsVxit) has an almost sinusoidal form with the 
dominating first and substantially weaker second frequency 
harmonics. The average drift velocity for the chosen parame- 
ters is vq/vf ~ 0.661 which corresponds to the fundamental 
frequency © /i = / ~ 0.661vF/a,x — 0.66 — 3.3 THz for 
ax — 1 — 0.2 fim. At S* ;» 1 the time dependence of the current 
is strongly non-monochromatic with quite large second, third, 
fourth and even sixth frequency harmonics. The average drift 
velocity for the chosen parameters is vo/vp « 0.615 with 
the resulting fundamental frequency /i = / w 0.615vp/ax- 
This corresponds, for the same values of the grating period 
Ux — 1 — 0.2 /im, to the first harmonic /i ~ 0.6 — 3 
THz, second harmonic /2 ~ 1.2 — 6 THz, fourth harmonic 
/4 ~ 2.4 — 12 THz and so on. As seen from Figure Ub) even 
the ninth harmonic (fg is up to 27 THz) has the amplitude 
only one order of magnitude smaller than the first one. The 
emitted radiation is linearly polarised in the direction of the 
dc current (x-direction). 

So far we have discussed the traditional field-effect- 
transistor-type design with the top graphene layer having only 
one (gate) contact. Figure |2}l. Alternatively, the top layer may 
have two metallic contacts, "source 2" and "drain 2", which 
gives the opportunity to drive a dc current in the second layer 
too. Figure |2^. Then the system becomes symmetric, with the 
second (first) layer serving as the grating coupler for the first 
(second) layer and the opportunity to emit radiation with both 
X- and ^/-polarization independently controlled by the source- 
drain voltages Vsdi and Vsd2- Another possible embodiment 
employing a single graphene layer is shown in Figure |2f . The 
periodic potential U{x) oc 1/Wy{x) is produced in this case 
by the alternating stripe width Wy{x). An advantage of this 
design is that it employs only one graphene layer In this case, 
however, the mean free path should exceed the grating period 

ax- 

In addition to the ability to work as an emitter, the proposed 
device may also operate as a new type of a field-effect- 
transistor (amplifier) combined with a plane radiating antenna. 
Indeed, the intensity of the emitted radiation is determined by 
the amplitude of the ac electric current ji. This amplitude 
depends, in its turn, on the voltage V12 between the first 
and second graphene layers: if 14 2 = 0, the current flowing 
from the source to the drain does not depend on time and no 
radiation is emitted; if V12 7^ 0, the current gets modulated 
and the device emits an electromagnetic signal. The radiation 
intensity can be very large (~ 5 kilowatts from square meter, 
Eq. ( [Tol l) at the few-millivolt input signal V12 (see estimates 
above). The proposed device thus amplifies an input signal 
and sends it directly to the surrounding space within the same 
physical process. 



The single-layer graphene absorbs only 2.3% of visible light 
ill 21 . The h-BN layer is a dielectric with a large band gap and 
is also transparent. The proposed powerful terahertz devices 
will therefore be very thin, light and practically invisible. 
Furthermore, such few-nm thick devices can be bent up 
focusing radiation and producing a huge concentration of THz 
power in a very small spatial volume. 

To summarise, we have proposed a new type of few- 
atomic-layers thin, light, bendable, almost invisible and volt- 
age tunable graphene based device which is able to produce 
a powerful electromagnetic radiation in a broad range of 
terahertz frequencies. This device may cause a revolution in 
the terahertz science and technology and will have a lot of 
applications in many different areas. 
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